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ABSTRACT 

We  present  an  analysis  of  X-ray  and  UV  data  obtained  with  the  XMM-Newton  Observatory  of  the  long-period  dwarf 
nova  RU  Peg.  RU  Peg  contains  a  massive  white  dwarf  (WD),  possibly  the  hottest  WD  in  a  dwarf  nova  (DN),  it  has  a 
low  inclination,  thus  optimally  exposing  its  X-ray  emitting  boundary  layer  (BL),  and  has  an  excellent  trigonometric 
parallax  distance.  We  modeled  the  X-ray  data  using  XSPEC  assuming  a  multi-temperature  plasma  emission  model 
built  from  the  MEKAL  code  (i.e.,  CEVMKL).  We  obtained  a  maximum  temperature  of  31.7  keV,  based  on  the 
European  Photon  Imaging  Camera  MOS 1 ,  2  and  pn  data,  indicating  that  RU  Peg  has  an  X-ray  spectrum  harder  than 
most  DNe,  except  U  Gem.  This  result  is  consistent  with  and  indirectly  confirms  the  large  mass  of  the  WD  in  RU  Peg. 

The  X-ray  luminosity  we  computed  corresponds  to  a  BL  luminosity  for  a  mass  accretion  rate  of  2  x  10-1 1  M0  yr_  1 
(assuming  Mw( j  =  1.3  M0),  in  agreement  with  the  expected  quiescent  accretion  rate.  The  modeling  of  the  O  vm 
emission  line  at  19  A  as  observed  by  the  Reflection  Grating  Spectrometer  implies  a  projected  stellar  rotational 
velocity  vrot  sin  i  —  695  km  s_1,  i.e.,  the  line  is  emitted  from  material  rotating  at  ~936-1245  km  s”1  (i  ~  34° -48°) 
or  about  1/6  of  the  Keplerian  speed;  this  velocity  is  much  larger  than  the  rotation  speed  of  the  WD  inferred  from 
the  Far  Ultraviolet  Spectroscopic  Explorer  spectrum.  Cross-correletion  analysis  yielded  an  undelayed  (time  lag 
~  0)  component  and  a  delayed  component  ofll6±17s  where  the  X-ray  variations/fluctuations  lagged  the  UV 
variations.  This  indicates  that  the  UV  fluctuations  in  the  inner  disk  are  propagated  into  the  X-ray  emitting  region  in 
about  116  s.  The  undelayed  component  may  be  related  to  irradiation  effects. 
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1.  INTRODUCTION 

1.1 .  The  Long-period  Dwarf  N ova  RU  Peg 

Dwarf  novae  (DNe)  are  a  class  of  weakly  magnetic  cata¬ 
clysmic  variables  (CVs)  which  are  interacting  compact  binaries 
in  which  a  white  dwarf  (WD,  the  primary  star)  accretes  mat¬ 
ter  and  angular  momentum  from  a  main-  (or  post-main-)  se¬ 
quence  star  (the  secondary)  filling  its  Roche  lobe.  The  matter 
is  transferred,  at  a  continuous  or  sporadic  accretion  rate  (M), 
by  means  of  an  accretion  disk  usually  reaching  all  the  way  to 
the  WD  surface.  Ongoing  accretion  at  a  low  rate  (quiescence)  is 
interrupted  every  few  weeks  to  months  by  intense  accretion  (out¬ 
burst)  of  days  to  weeks  (a  DN  accretion  event).  DNe  are  powerful 
X-ray  sources  with  luminosities  of  1030— 1033  erg  s  .  The 
X-ray  emission  is  believed  to  originate  in  the  boundary  layer 
(BL)  between  the  slowly  rotating  accreting  WD  and  the  fast 
rotating  (Keplerian)  inner  edge  of  the  accretion  disk,  where  the 
material  dissipates  its  remaining  rotational  kinetic  energy  before 
it  accretes  onto  the  surface  of  the  WD.  The  typical  DNe,  i.e., 
those  systems  exhibiting  normal  DN  outbursts,  are  the  U  Gem 
sub-type  of  DNe  (according  to  the  classification  of  Ritter  & 
Kolb  2003)  and  are  located  above  the  period  gap  (Porb  >  3  hr). 

RU  Peg  is  a  U  Gem  type  DN  with  an  orbital  period  Porb  = 
8.99  hr,  a  secondary  spectral  type  K2-5V,  a  primary  (WD) 
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mass  Mwd  =  1.29*  qjq  Mq,  and  a  secondary  mass  M2  — 
0.94  ±  0.04  M0  (Stover  1981;  Wade  1982;  Shafter  1983).  The 
system  has  a  magnitude  range  Umax- Emin  ~  9.0-13.1  with 
outbursts  lasting  ~20  days  and  recurring  every  ~50  days.  The 
near-Chandrasekhar  WD  mass  has  been  corroborated  by  the 
sodium  (8190  A)  doublet  radial  velocity  study  of  Friend  et  al. 
(1990).  They  obtained  a  mass  of  1.38  ±  0.06  M0  for  the  WD  and 
also  found  a  range  of  inclination  angles  between  34°  and  48°  in 
agreement  with  the  range  of  plausible  inclinations  found  in  the 
study  by  Stover  (1981).  More  recently,  a  Hubble  Fine  Guidance 
Sensor  (FGS)  parallax  of  3.55  ±  0.26  mas  was  measured  by 
Johnson  et  al.  (2003)  yielding  a  distance  of  282  ±  20  pc. 

RU  Peg  was  observed  with  IUE  under  several  different 
observing  programs,  both  in  quiescence  and  in  outburst,  and 
was  part  of  several  survey-like  studies  (e.g.,  La  Dous  et  al. 
1985;  Verbunt  1987;  Szkody  et  al.  1991  to  cite  just  a  few). 

Sion  &  Urban  (2002)  modeled  four  IUE  spectra  obtained  in 
deep  quiescence  with  accretion  disks  and  photospheres.  They 
found  that  a  very  hot  WD  dominated  the  FUV  spectrum  with  a 
temperature  Peff  =  50-53,000  K  which  places  RU  Peg  among 
the  hottest  WDs  in  DNe.  The  distance  corresponding  to  their 
best  fitting,  high  gravity  (Log(g)  =  8.7)  photosphere  models 
was  250  pc. 

More  recently,  Godon  et  al.  (2008)  modeled  the  Far  Ultra¬ 
violet  Spectroscopic  Explorer  (FUSE)  spectrum  of  RU  Peg  in 
quiescence  and  obtained  a  WD  with  a  temperature  of  70,000  K, 
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a  rotational  velocity  of  40  km  s_1,  assuming  Log(g)  =  8.8  and 
a  distance  of  282  pc.  In  this  later  study,  the  higher  tempera¬ 
ture  obtained  in  the  model  fitting  is  mainly  a  consequence  of 
the  assumed  larger  distance  and  gravity.  It  is  clear  that  RU  Peg 
has  a  massive  WD  and  therefore  a  deep  potential  well,  and  its 
surface  temperature  is  very  large  (>50,000  K),  possibly  point¬ 
ing  to  strong  accretional  and  BL  heating.  For  these  reasons,  we 
chose  RU  Peg  as  our  X-ray  target,  as  it  is  expected  to  be  a  copi¬ 
ous  source  of  X-rays  and  should  be  an  ideal  candidate  to  study 
its  BL. 

1.2.  The  Boundary  Layer 

The  BL  is  the  region  between  the  slowly  rotating  accreting 
WD  and  the  fast  rotating  (Keplerian)  inner  edge  of  the  accretion 
disk.  For  accretion  to  occur,  gravity  has  to  overcome  the  cen¬ 
trifugal  acceleration,  and  this  happens  in  the  BL  as  the  material 
dissipates  its  remaining  rotational  kinetic  energy  before  being 
accreted  onto  the  surface  of  the  WD.  In  the  BL  the  rotational 
velocity  is  sub-Keplerian  and  decreases  inward. 

The  standard  disk  theory  (Shakura  &  Sunyaev  1973)  predicts 
a  BL  luminosity  almost  equal  to  the  disk  luminosity  (Pringle 
1981): 

Lbl  =  {1- P2)Ldisk  =  (l- p2)—^,  (1) 

where  G  is  the  gravitational  constant,  M*  is  the  mass  of 
the  accreting  star,  R,,  is  its  radius,  M  is  the  mass  accretion 
rate,  and  ft  is  the  stellar  angular  velocity  in  Keplerian  units 
P  —  7?*).  For  most  systems  the  WD  rotational  velocity 

is  of  the  order  of  a  few  100  km  s-1  and  one  therefore  has 
P  <<£  1  and  Lbl  ~  7-disk-  For  a  star  rotating  at  (e.g.)  10% 
of  the  breakup  velocity  Equation  (1)  gives  Lbl  =  0. 997-disk- 
However,  the  standard  disk  theory  does  not  take  into  account 
the  BL  explicitly.  Using  a  one-dimensional  approach,  Kluzniak 
(1987)  has  shown  that  part  of  the  BL  (kinetic)  energy  actually 
goes  into  the  spinning  of  the  star  in  the  equatorial  region,  and 
the  BL  luminosity  is  more  accurately  given  by  the  relation 

7-bl  =  (1  —  P)1  7-disk-  (2) 

For  a  star  rotating  at  10%  of  the  breakup  velocity,  Equation  (2) 
gives  7-bl  =  0.81Ldisk,  which  is  significantly  different  than 
Equation  (1). 

Because  of  its  small  radial  extent,  the  BL  is  expected  to  emit 
its  energy  in  the  X-ray  bands  (7-bl  ^  7-x-ray)-  At  high  accretion 
rate  (or  during  DN  outburst),  the  BL  is  expected  to  be  optically 
thick  with  a  temperature  T  ~  105— 106  K  (Godon  et  al.  1995; 
Popham  &  Narayan  1995)  and  emits  in  the  soft  X-ray  band. 
Observations  of  CVs  in  high  state  (e.g.,  Mauche  et  al.  1995; 
Baskill  et  al.  2005)  confirm  these  predictions.  At  low  mass 
accretion  rate  (or  during  DN  quiescence),  as  the  density  is  much 
decreased,  the  BL  becomes  optically  thin  and  emits  in  the  hard 
X-ray  band  with  T  &  108  K  (Narayan  &  Popham  1993)  (see 
below  for  observational  evidence  of  such  optically  thin  BLs). 
Hence,  during  quiescence  the  emission  should  arise  from  a  very 
hot  plasma  very  close  to  the  WD  surface. 

Previous  X-ray  observations  of  DN  systems  in  quiescence 
(e.g.,  van  der  Woerd  &  Heise  1987;  van  Teeseling  &  Verbunt 
1996;  Belloni  et  al.  1991),  while  confirming  the  presence  of 
hard  X-ray,  deduced  that,  contrary  to  the  theory,  the  quies¬ 
cent  BL  luminosity  was  underluminous.  Namely,  they  found 
LBl  =  7-x-ray  7-disk,  which  confirmed  the  original  claim 


that  BLs  were  actually  missing  (Ferland  et  al.  1982).  How¬ 
ever,  these  earlier  results  assumed  the  disk  to  be  the  source 
of  the  optical  and  ultraviolet  radiation  (Ldisk  =  Lopt  +  Luv) 
and  used  Equation  (1)  rather  than  Equation  (2).  More  recently, 
X-ray  Multi  Mirror-Newton  ( XMM-Newton )  observations  of 
eight  DNe  in  quiescence  by  Pandel  et  al.  (2003,  2005)  re¬ 
vealed  that  a  significant  part  of  the  emitted  FUV  flux  (Luv) 
actually  originates  from  the  WD  itself,  and  the  evidence  for  un¬ 
derluminous  BLs  in  quiescent  DNe  was  refuted.  For  RU  Peg, 
FUSE  observations  (Godon  et  al.  2008)  indicate  that  the  WD 
contributes  possibly  most  of  the  FUV  light  with  a  temperature 
T  >  50,000  K.  In  addition,  Godon  &  Sion  (2005)  noted  that 
the  region  where  the  BL  meets  the  inner  edge  of  the  Keplerian 
disk  (Popham  1999)  can  also  contribute  some  FUV  flux.  It  is 
clear  now  that  one  cannot  just  compare  the  X-ray  luminosity 
to  the  optical  +  UV  luminosity  to  check  whether  the  boundary 
luminosity  is  as  large  as  the  disk  luminosity. 

For  the  eight  DNe  caught  in  quiescence,  Pandel  et  al. 
(2003,  2005),  using  the  XMM-Newton  data,  obtained  X-ray 
BL  luminosities  of  the  order  of  ~1  x  1031  erg  s'1  to  6.6  x 
1032  erg  s-1,  with  temperatures  ranging  from  ~8  to  55  keV,  and 
mass  accretion  rates  deduced  from  X-rays  in  the  range  1 01 2  Me 
yr”1  to  10~10  Me  yr”1.  In  this  work,  we  present  XMM-Newton 
observations  of  RU  Peg  taken  in  quiescence  to  derive  its  X-ray 
luminosity  and  gain  information  on  its  BL. 

2.  OBSERVATIONS  AND  ANALYSIS 

Th e.  XMM-Newton  Observatory  (Jansen  et  al.  2001)  has  three 
1500  cm2  X-ray  telescopes  each  with  an  European  Photon 
Imaging  Camera  (EPIC)  at  the  focus.  Two  of  the  telescopes 
have  Multi-Object  Spectrometer  (MOS)  CCDs  (Turner  et  al. 
2001)  and  the  last  one  uses  pn  CCDs  (Striider  et  al.  2001)  for 
data  recording.  There  are  two  Reflection  Grating  Spectrometers 
(RGSs;  den  Herder  et  al.  2001).  The  Optical  Monitor  (OM),  a 
photon  counting  instrument,  is  a  co-aligned  30  cm  optical/UV 
telescope,  providing  for  the  first  time  the  possibility  to  observe 
simultaneously  in  the  X-ray  and  optical/UV  regime  from  a 
single  platform  (Mason  et  al.  2001). 

RU  Peg  was  observed  (pointed  observation)  with  XMM- 
Newton  for  a  duration  of  53.1  ks  on  2008  June  9  at  07:16:50.0 
UTC  (obsID  0551920101).  At  that  time  the  system  was  at  a 
visual  magnitude  of  ~  1 2.5,  about  2  months  into  quiescence 
and  2  weeks  before  the  next  outburst  (from  AAVSO  data8). 
Data  were  collected  with  the  EPIC  MOS  and  pn  cameras  in  the 
prime  partial  window2  and  prime  full  window  imaging  mode, 
respectively,  the  RGS  and  the  OM  using  the  fast  imaging  mode 
(L0.5  s  time  resolution)  with  the  UVW1  filter  (240-340  nm). 

We  analyzed  the  pipeline-processed  data  using  Science  Anal¬ 
ysis  Software  (SAS)  version  9.0.0.  Data  (single-  and  double¬ 
pixel  events,  i.e.,  patterns  0-4  with  Flag  =  0  option  for  pn  and 
patterns  ^12  with  Flag  =  0  for  MOS  1,2)  were  extracted  from 
a  circular  region  of  radius  60"  for  pn  and  45"  MOS  1  and  MOS2 
in  order  to  perform  spectral  analysis  together  with  the  back¬ 
ground  events  extracted  from  a  source  free  zone  normalized  to 
the  source  extraction  area.  We  checked  the  pipeline-processed 
event  file  for  any  existing  flaring  episodes  and  no  sporadic  events 
in  the  background  were  detected  with  count  rate  higher  than 
0.08  counts  s_1  (for  MOS  1,2)  and  0.5  for  pn  detectors.  Table  1 
displays  the  background  subtracted  count  rates  for  the  EPIC  pn, 
MOS1,  and  MOS2. 
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Time(s) 

Figure  1.  EPIC  X-ray  light  curve  (binned  at  16  s,  in  red)  together  with  the  OM  UV  light  curve  (binned  here  at  100  s  for  display,  in  black;  the  data  were  binned  at  20  s 
for  the  correlation  check).  The  count  rate  of  the  UV  data  has  been  divided  by  9. 1 736  to  fit  the  count  rate  level  of  the  X-ray  data  for  easier  comparison.  The  time  t  =  0 
corresponds  to  t  =  3.29389320  x  108  s,  counted  from  the  MJD  reference  day  50,814,  namely  54,626.376  MJD.  The  time  modulation  of  the  UV  data  follows  closely 
the  time  modulation  of  the  X-ray  data  except  around  (  at  12  ks,  20-22  ks,  31  ks,  and  48.5  ks,  where  the  UV  has  relatively  more  flux. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Table  1 

XMM-Newton  Observations  of  RU  Peg 


Instrument 

Count  Rate 

(s-1) 

RGS2 

0.2881  ±  0.0030 

RGS1 

0.2351  ±  0.0029 

EPIC  pn 

10.03  ±0.017 

EPIC  MOS 1 

3.154  ±0.0084 

EPIC  MOS2 

3.241  ±  0.0085 

The  RGS  observations  were  carried  out  using  the  standard 
spectroscopy  mode  for  readout.  We  reprocessed  the  data  using 
the  XMM-SAS  routine  rgsproc.  We  first  made  event  files 
and  determined  times  of  low  background  from  the  count  rate 
on  CCD  9  (which  is  closest  to  the  optical  axis).  The  final 
exposure  times  and  net  count  rates  showed  that  there  were  no 
sporadic  high  background  events  in  our  data.  Table  1  displays  the 
background-subtracted  count  rates  for  RGS  1  and  RGS2.  Source 
and  background  counts  for  the  RGS  were  extracted  using  the 
standard  spatial  and  energy  filters  for  the  source  position,  which 


defines  the  spatial  extraction  regions  as  well  as  the  wavelength 
zero  point. 

2.1.  X-Ray  and  UV  Light  Curves 

The  UV  and  EPIC  pn  X-ray  light  curves  are  shown  together 
for  comparison  in  Figure  1,  where  we  have  scaled  the  UV  count 
rate  to  fit  the  X-ray  count  rate  for  a  better  comparison.  The 
count  rate  for  the  UV  data  ranges  between  ~50  counts  s-1  and 
~100  counts  s~  1  with  a  time  average  of  ~7 1 .8  counts  s-1.  For 
UV-bright  objects,  a  count  rate  of  1  s~  1  in  the  UVW1  filter 
translates  into  a  flux  of  4.5  x  10~16  erg  cm~2  s-1  A-1  at  290 
nm  (see  the  online  XMM-Newton  documentation9 *).  This  gives 
for  RU  Peg  a  flux  of  ^2. 2-4. 5  x  10~14  erg  cm-2  s_I  A-1, 
corresponding  to  a  luminosity  of  ^(2. 1-4.3)  x  1029  erg  s-1  A'1 
(at  a  distance  d  =  282  pc).  The  time  modulation  of  the  UV  data 
follows  closely  the  time  modulation  of  the  X-ray  data  except 
around  t  12  ks,  20-22  ks,  31  ks,  and  48.5  ks,  where  the 
UV  has  relatively  more  flux  for  a  duration  of  several  hundred 
seconds  (and  up  to  1000  s).  Since  the  UV  is  expected  to  be 


9  http://xmm.vilspa.esa.es/extemal/xmm_user_support/documentation/ 

uhb/index.html 
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Time  Delay  (s)  Time  De|aY  (s) 

Figure  2.  Cross-correlation  between  X-ray  and  UV  light  curves.  Left:  a  single  Lorentzian  is  used  to  fit  the  data  at  0;  right:  a  second  Lorentzian  fit  is  added  to  fit  the 
delayed  component.  The  figures  show  strong  correlation  at  zero  time  lag.  The  asymmetric  profile  centered  at  zero  time  lag  (asymmetry  toward  positive  lag)  indicates 
that  the  X-ray  variations  are  delayed  with  respect  to  the  UV  variations.  Error  bars  indicate  the  standard  deviations  from  the  average  of  the  value  of  the  cross-correlation 
at  any  given  lag  in  different  time  segments  where  power  spectra  are  calculated. 


emitted  further  out  than  the  X-rays,  these  four  epochs  where  the 
UV  light  curves  do  not  decrease  as  much  as  the  X-rays  might  be 
due  to  the  occultation  of  the  X-ray  emitting  material  by  the  WD, 
while  the  UV  emitting  region  is  not  hidden  from  the  observer. 

In  order  to  study  the  correlation  between  the  X-ray  and  the  UV 
variability,  we  calculated  the  cross-correlation  between  the  two 
light  curves.  We  used  time  bins  of  5  s  averaging  several  power 
spectra  with  128  bins  for  the  analysis.  The  resulting  correlation 
coefficient  as  a  function  of  time  lag  is  shown  in  Figures  2(a)  and 
(b).  The  correlation  coefficient  is  normalized  to  have  a  maximum 
value  of  1.  The  curve  shows  a  clear  asymmetry  indicating  the 
existence  of  time  delays.  We  also  detect  a  strong  peak  near  zero 
time  lag  suggesting  a  significant  correlation  between  X-rays  and 
the  UV  light  curves.  We  call  this  the  undelayed  component  (see 
Figures  2(a)  and  (b)).  The  positive  time  lag  in  the  asymmetric 
profile  shows  that  the  X-ray  variations  are  delayed  relative  to 
those  in  the  UV.  In  order  to  calculate  an  average  time  lag  that 
would  produce  the  asymmetric  profile,  we  fitted  the  varying 
cross-correlation  by  two  Lorentzians,  with  time  parameter  fixed 
at  0.0  lag  and  the  other  set  as  free.  The  resulting  fit  yields  a  lag 
ofll6±  17  s.  This  is  the  delayed  component. 

2.2.  EPIC  Spectrum 

We  performed  spectral  analysis  of  the  EPIC  data  using  the 
SAS  task  ESPECGET  and  derived  the  spectra  of  the  source 
and  the  background  together  with  the  appropriate  response 
matrices  and  ancillary  files.  How  the  photons  were  extracted 
is  described  in  Section  2.  The  EPIC  pn,  EPIC  MOS1,  and 
EPIC  MOS2  spectra  were  simultaneously  fitted  to  derive  the 
spectral  parameters.  The  spectral  analysis  was  performed  using 
XSPEC  version  12.6.0q  (Arnaud  1996).  A  constant  factor  was 
included  in  the  spectral  fitting  to  allow  for  a  normalization 
uncertainty  between  the  EPIC  pn  and  EPIC  MOS  instruments. 
We  grouped  the  pn  and  MOS  spectral  energy  channels  so  that 
there  is  a  minimum  of  80  (MOS  1,2)  to  150  (pn)  counts  in  a 
bin  to  improve  the  statistical  quality  of  the  spectra.  The  fits 
were  conducted  in  the  0.2-10.0  keV  range.  The  simultaneously 
fitted  spectra  from  the  three  EPIC  instruments  are  shown  in 
Figure  3.  We  modeled  the  X-ray  spectrum  of  RU  Peg  in  a 
similar  fashion  as  Pandel  et  al.  (2005)  and  fitted  the  data 
with  the  (TBabs  x  CEVMKL)  model  within  XSPEC.  TBabs 
is  the  Tuebingen-Boulder  ISM  absorption  model  (Wilms  et  al. 
2000)  and  CEVMKL  is  a  multi-temperature  plasma  emission 
model  built  from  the  mekal  code  (Mewe  et  al.  1985).  Emission 


Figure  3.  EPIC  MOS  1,2  and  pn  spectra  simultaneously  fitted  to  the  same 
model  TBabs  x  CEVMKL.  The  model  includes  an  ISM  absorption  model  and 
a  multi-temperature  plasma  emission  model  built  from  the  MEKAL  code.  The 
maximum  shock  temperature  obtained  from  the  fit  is  3 1 .7  keV.  The  second  panel 
below  the  fitted  spectra  is  the  residuals  in  standard  deviations.  One  can  notice  a 
slight  CTI  problem  of  the  pn  detector  around  the  iron  lines  (6.4-6. 9  keV)  where 
the  fluctuations  in  the  residuals  due  to  CTI  are  evident.  This  does  not  affect  the 
global  fitting  procedure  but  only  slightly  increases  the  reduced  x2  of  the  fits. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

measures  follow  a  power  law  in  temperature  (i.e.,  emission 
measure  from  temperature  T  is  proportional  to  (T /  7max)"  1 ). 
The  residuals  in  Figure  3  show  systematic  fluctuation  around 
the  6. 7-6. 9  keV  iron  line  complex  mainly  from  the  EPIC  pn 
data  and  some  small  low-energy  fluctuations  exist  in  the  MOS2 
data  as  well.  This  is  due  to  the  charge  transfer  inefficiency 
(CTI)  problem  in  the  pn  (and  possibly  MOS2)  instrument. 
These  generally  occur  around  lines  due  to  small  calibration 
errors  and  mostly  effect  only  the  line  shapes  leaving  systematic 
residuals  and  increasing  the  reduced  x2  of  the  fits.  Our  EPIC 
MOS1  data  do  not  exhibit  any  CTI  effects  and  the  reduced 
X2  for  the  fit  to  these  data  alone  is  1.15  (dof  290).  The 
reduced  x 2  for  the  simultaneously  fitted  spectra  is  higher  than 
the  value  for  the  MOS1  fit,  but  the  spectral  parameters  for 
all  three  instruments  are  almost  the  same  within  the  errors. 
Table  2  contains  the  spectral  parameters  from  fits  (using  three 
detectors  simultaneously)  with  the  (TBabs  x  CEVMKL)  model. 
Errors  are  given  at  the  90%  confidence  level.  We  find  a 
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Table  2 

Spectral  Parameters  of  Fit  to  the  Combined  EPIC  Spectrum  of  RU  Peg  in  the 
Energy  Range  0.2-10  keV 


Parameter 

Value 

Nn  (1022  atoms  cm'2) 

0  044+0  001 
u.u^_00oi 

a 

1  05+003 
i-UD-0.03 

Tmax  (keV) 

31  7+17 
di.  /_20 

o 

ft  72+0.06 
u'  — 0.06 

Ne 

n  55+O.I6 
u-:o—  0.16 

Mg 

1  3+0‘2 

A,:,-0.2 

Si 

ft  o+0.14 

0.5_o  i4 

S 

n  q+0.3 

U-  — 0.2 

Ca 

i  q+0.8 

0.7 

Fe 

r»  o+0.04 

0-O— 0.04 

^CEVMKL 

0  o47+0  001 
u.U‘t/_oooi 

Gaussian  LineE  (keV) 

6.4  (fixed) 

aG  (keV) 

a  i  q+0.025 

U,A  ^-0.025 

Kg 

o  oooo59+0  000006 

u.uuuudz_0  000006 

Flux  (ten11  erg  cm-2  s-1) 

4  1+0.2 

Xv  Wot) 

1.5  (1741) 

Notes.  Nn  is  the  absorbing  column,  a 

is  the  index  of  the  power-law 

emissivity  function  (dEM  =  ( T/Tmax)a  ^dT/Tm axX  7max  is  the 
maximum  temperature  for  the  CEVMKL  model.  Element  names 

stand  for  the  abundance  relative  to 

solar  abundances,  Gaussian 

lineE  is  the  line  center  for  the  emission  line,  cxq  is  the  line  width; 
^cevmkl  and  Kq  are  the  normalizations  for  CEVMKL  and  Gaussian 
models,  respectively.  The  unabsorbed  X-ray  flux  is  given  in  the  range 
0.2-10.0  keV.  All  error  ranges  are  given  at  the  90%  confidence  level 

(Ax2  =  2.71  for  a  single  parameter). 

Table  3 

Line  Identifications  of  the  XMM-Newton  RGS  Spectrum  of  RU  Peg 

Ion 

Wavelength 

(A) 

Mg  xn 

7.20 

7.45 

Nex 

12.15 

Fe  xvn 

15.03 

O  viii 

16.05 

Fe  xvii 

16.80 

Fe  xvii 

17.07 

O  viii 

19.00 

O  vn  (i) 

21.80 

O  vn  (f) 

22.11 

N  vii 

24.77 

maximum  plasma  temperature  in  a  range  29-33  keV  and 
mostly  solar  abundances  of  elements  aside  from  oxygen  and 
neon  which  we  calculate  to  be  subsolar.  The  unabsorbed 
X-ray  flux  is  4.1^j22  x  10~n  erg  s^1  cirT2  which  translates  to  a 
luminosity  of  4.Lq3  x  1032  erg  s-1  at  282  pc  (see  Section  1). 
The  neutral  hydrogen  column  density  is  (4. 3-4. 5)  x  1020  cm~2. 

2.3.  RGS  Spectrum 

The  RGS  analysis  tool  rgsproc  was  used  to  obtain  RGS1 
and  RGS2  spectra  and  to  produce  a  fluxed  spectrum  (i.e.,  using 
rgsfluxer).  The  resultant  fluxed  spectrum  of  the  combined 
RGS  1  and  RGS2  detectors  is  shown  in  Figure  4  with  line  identi¬ 
fications.  The  detected  lines  and  corresponding  wavelengths  are 
listed  in  Table  3.  Fluxed  spectra  are  obtained  just  by  dividing  the 


1.4 


30  31  32  33  34  35 

X  (A) 

Figure  4.  Fluxed  RGS  Spectrum  of  RU  Peg  with  line  identifications. 


count  spectrum  by  the  RGS  effective  area.  It  neglects  the  redistri¬ 
bution  of  monochromatic  response  into  the  dispersion  channels. 
Since  proper  response  is  not  utilized  by  the  fluxed  spectrum,  in 
order  to  perform  a  similar  spectral  fit/analysis  using  the  same 
(TBabs  x  CEVMKL)  model  within  XSPEC,  we  used  the  count 
rate  spectra  produced  for  RGS1  and  RGS2  simultaneously  with 
the  appropriate  response  files  for  each  detector.  This  is  a  very 
efficient  approach  to  find  the  spectral  parameters  in  comparison 
with  the  EPIC  results.  The  fitted  RGS1  and  RGS2  spectra  are 
shown  in  Figure  5.  Table  4  contains  the  spectral  parameters  from 
the  fit  with  the  (TBabs  x  CEVMKL)  model  within  XSPEC. 
Errors  are  at  the  90%  confidence  level  and  the  fit  is  performed 
between  0.2  and  2.5  keV. 

The  maximum  temperature  from  the  fits  with  the  RGS  data 
yields  24  keV  with  a  large  error  range  of  17 — 41  keV  since 
the  spectrum  has  a  lower  count  rate  and  a  lower  upper  energy 
boundary  (i.e.,  2.5  keV)  compared  with  the  EPIC  data  (i.e., 
12  keV).  We  find  that  most  of  the  RGS  spectral  parameters  are 
consistent  with  the  EPIC  results.  Oxygen  and  neon  are  subsolar 
in  abundance  and  silicon  appears  to  be  slightly  enhanced 
compared  to  solar  abundance.  The  resolution  of  the  RGS 
spectra  is  sufficient  to  measure  the  rotational  velocity  of  the 
BL  via  Doppler  broadening  of  emission  lines.  We  calculated 
the  broadening  using  the  O  vm  emission  line  at  19  A  which  is 
the  strongest  line  in  the  spectrum.  We  used  a  Gaussian  model 
to  calculate  the  a  of  the  line  (a  x  2.4  =  FWHM)  along  with 
a  power  law/bremmstrahlung  for  the  continuum.  We  find  that 
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Figure  5.  Simultaneously  fitted  RGS1  and  RGS2  spectrum  of  RU  Peg.  The 
same  (TBabsxCEVMKL)  model  is  used  for  the  fit.  The  residuals  in  the  second 
panel  below  are  in  standard  deviations. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

Table  4 

Spectral  Parameters  of  the  Fit  to  the  RGS  Spectrum  of  RU  Peg  in  the  Energy 
Range  0.2-2. 5  keV  with  Same  Definitions  as  in  Table  2 


Parameter 

Value 

Nn  (1022  atoms  crrT2) 

0  O31+0  007 
U.Uj1_q  007 

a 

1  9+O.I6 
i>z-0.16 

Lnax  (keV) 

94  1+17.0 
z^'A-7.0 

O 

r»  O+0.13 
U.O—o  12 

Ne 

I  1+0.5 

1 '  — 0.5 

Mg 

2  7+1-2 
z-'— 1.1 

Si 

o  4+2.0 

3 -^-2.0 

S 

1  4< 

1  •  — 3.2 

Ca 

1 .0  (fixed) 

Fe 

1.0  (fixed) 

S'CEVMKL 

0  043+0  004 

u.uho_q  005 

Flux  (10~n  erg  s-1  cm-2) 

t  q+0.4 
~’-y-0.4 

Xv2  (dof) 

1.3  (283) 

FWHM  =  0.044  A  and  the  corresponding  velocity  in  the  line 
of  sight  is  695  km  s-1  (A/.//.  =  v/c).  We  also  checked  the 
resolution  of  RGS  at  around  19  A  and  found  0.04-0.05  A  and 
we  caution  that  our  measurement  is  on  the  limit  of  the  RGS 
spectral  resolution. 


3.  DISCUSSION 

The  maximum  X-ray  shock  temperature  we  obtain  is 
29-33  keV,  based  on  the  EPIC  MOSl,2  and  pn  data,  as  they 
have  higher  count  rates  and  broader  energy  ranges  than  the 
RGS  data.  RU  Peg  has  a  harder  X-ray  spectrum  than  in  most 
DNe,  but  softer  than  U  Gem  in  quiescence,  which  also  contains 
a  massive  WD  (Sion  et  al.  1998;  Long  &  Gilliland  1999).  This 
result  is  consistent  with  and  confirms  the  large  mass  of  the  WD 
in  RU  Peg. 

The  X-ray  luminosity  of  RU  Peg  is  4.1  x  1032  erg  s~  1 
(0.2-10.0  keV),  and  assuming  we  see  only  half  of  the  BL  (if 
it  is  close  to  the  star),  we  have  Lbl  =  8.2  x  1032  erg  s  '. 
This  is  in  the  range  of  the  other  quiescent  DNe  observed  with 
XMM-Newton  by  Pandel  et  al.  (2005).  However,  in  order  to 


Figure  6.  Boundary  layer  (X-ray)  luminosities  vs.  the  UV  spectral  luminosity 
at  290  nm,  as  derived  by  Pandel  et  al.  (2005)  for  quiescent  DNe  observed  with 
XMM-Newton.  We  have  added  RU  Peg  on  the  graph.  The  solid  and  dotted 
lines  show  Lbl  =  Em  (=  ^uvwi)  assuming  UV  luminosity  predicted  by  a 
simple  accretion  disk  model  for  an  inner  radius  Rm  =  5000  km  and  10,000  km, 
respectively.  For  systems  higher  above  the  line,  such  as  RU  Peg,  the  UV  excess 
is  due  to  the  contribution  from  the  WD.  The  UV  luminosity  of  a  WD  with  an 
8000  km  radius  at  various  temperature  is  shown  on  the  left.  RU  Peg  corresponds 
here  to  a  ~53,350  K  WD  with  a  radius  of  8000  km  (which  we  have  also  added  to 
the  original  graph),  which  translates  to  ~75,000  K  for  a  4000  km  WD,  more  in 
line  with  the  large  mass  of  the  WD  in  RU  Peg.  WX  Hyi  was  caught  in  outburst 
explaining  its  UV  excess. 


fully  compared  RU  Peg  with  these  other  systems,  we  also 
need  to  consider  the  UV  luminosity.  For  the  UV  we  use  the 
spectral  luminosity  at  290  nm  Luvwi  obtained  from  the  OM 
data  (see  Section  2.1).  We  have  reproduced  Figure  4  of  Pandel 
et  al.  (2005)  in  Figure  6  with  the  inclusion  of  RU  Peg,  which 
shows  the  quiescent  DNe  observed  with  XMM-Newton  plotted 
on  an  Luvwi  against  Lbl  graph.  In  this  figure  the  solid  and 
dotted  lines  show  the  relation  Lbl  =  Ldisk  for  a  UV  luminosity 
predicted  by  a  simple  accretion  disk  model  with  an  inner  radius 
of  5000  Km  and  10,000  Km,  respectively.  This  simple  disk 
model  does  not  include  UV  contribution  from  the  WD  or  outer 
edge  of  the  BL  where  it  meets  the  inner  disk.  The  location  of 
a  system  in  the  vicinity  of  this  diagonal  (e.g.,  SU  UMa,  WW 
Hyi)  indicates  that  the  BL  luminosity  Lbl  is  comparable  to  the 
disk  luminosity  (here  L^sk  =  Luvwi)-  For  RU  Peg,  as  for  U 
Gem  for  example,  the  excess  of  UV  is  due  to  the  contribution 
from  the  WD  marked  on  the  left  of  the  graph.  RU  Peg  has  a  UV 
luminosity  corresponding  to  53,350  K  for  an  8000  km  radius 
WD  (assuming  the  flux  scales  simply  as  oc  L4)  or  75,450  K  WD 
for  a  4000  km  radius  (more  consistent  with  the  large  mass  of 
RU  Peg).  This  confirms  the  temperature  of  the  WD  as  derived 
from  the  FUSE  spectrum  (Godon  et  al.  2008)  and  puts  RU  Peg 
(literally)  in  line  with  all  the  other  quiescent  DNe  such  as  VW 
Hyi,  U  Gem,  SU  UMa,  OY  Car,  and  AB  Dra,  on  the  Luvwi 
versus  Lbl  graph. 

Because  of  its  large  mass  (and  therefore  small  radius),  the  WD 
in  RU  Peg  has  a  deeper  potential  well,  and  i>Kep  of  the  order  of 
5-6000  km  s~  1  in  the  BL.  As  the  matter  is  decelerated  in  the  BL, 
the  X-ray’s  emitting  gas  has  velocities  of  a  few  1000  km  s_1. 
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The  modeling  of  the  Ovm  emission  line  at  19  A  implies  a 
projected  rotational  broadening  vTOt  sin/  =  695  km  s'1,  i.e.,  the 
line  is  emitted  from  material  rotating  at  ~  936-1245  km  s'1 
(since  i  ~  34° -48°)  or  about  1/6  of  the  Keplerian  speed.  This 
velocity  is  still  much  larger  than  the  rotation  speed  of  the  WD 
inferred  from  the  FUSE  spectrum  (40  km  s'1;  Godon  et  al. 
2008).  This  implies  that  the  X-ray  emission  comes  directly 
from  the  decelerating  BL  material.  It  is  possible  that  the  X-ray 
emission  originates  in  the  equatorial  region  of  the  WD  as  shown 
by  Piro  &  Bildsten  (2004),  who  found  that  poloidal  motion  may 
be  negligible  at  low  mass  accretion  rates  characteristic  of  DNe 
in  quiescence.  In  such  a  case  most  of  the  dissipated  energy 
is  radiated  back  into  the  disk,  which  may  justify  the  use  of  a 
one-dimensional  treatment  (such  as  Narayan  &  Popham  1993; 
Popham  1999). 

The  X-ray  luminosity  we  computed  corresponds  to  a  BL 
luminosity  (Equation  (2))  for  a  mass  accretion  rate  of  2.0  x 
10'11  Mq  yr'1  assuming  a  1.29  Me  WD  mass,  and  it  increases 
to  3.0  x  1  O' 1 1  Mq  yr'1  for  a  1.2  MQ  WD  mass.  This  is  entirely 
consistent  with  a  quiescent  accretion  rate. 

We  checked  the  correlation  between  the  variability  of  the 
X-ray  data  and  the  UV  data  and  found  two  components:  one 
delayed  and  the  other  undelayed.  The  X-ray  and  UV  emission 
originate  from  distinct  regions  in  the  binary.  Since  RU  Peg  is 
a  non-magnetic  system,  the  X-ray  emission  is  entirely  from 
the  inner  edge  of  the  BL  very  close  to  the  WD  surface.  As 
to  the  UV  radiation,  it  is  emitted  by  the  heated  WD,  the  very 
inner  disk,  and  also  that  region  where  the  outer  BL  meets  the 
disk  (Popham  1999).  The  significant  modulation  correlation  at 
At  ~  0  lag  is  expected  to  be  caused  by  reprocessing  of  X-rays 
(i.e.,  irradiation  by  X-rays)  in  the  accretion  disk.  Such  time  lags 
are  on  the  order  of  milliseconds  and  proportional  to  light  travel 
time  which  is  well  beyond  the  time  resolution  in  our  light  curves. 
The  delayed  component  of  At  ~  116  s  lag  is  much  longer  and 
cannot  be  produced  by  light  travel  effects  nor  by  reprocessing  of 
the  X-ray,  since  the  X-ray  trails  behind  the  UV.  The  only  viable 
explanation  is  that  the  time  lag  is  the  time  it  takes  for  matter 
to  move  inward  from  the  very  inner  disk  (emitting  in  the  UV) 
onto  the  stellar  surface  (emitting  in  the  X-ray).  This  is  the  time 
it  takes  to  spin  down  the  material  in  the  BL,  rspin  =  At  —  116  s. 
The  modulations  of  the  UV  component  and  its  lagging  X-ray 
counterpart  are  due  to  modulations  in  M.  A  comparable  time- 
delayed  (At  ~  100  s)  component  was  also  detected  for  VW  Hyi 
(Pandel  et  al.  2003)  and  a  much  shorter  one  (~  7  s;  Revnivtsev 
et  al.  2011)  was  detected  for  the  intermediate  polar  EX  Hya 
(indicating  that  the  transit  of  matter  through  magnetic  field  lines 
to  the  poles  is  faster  than  through  the  non-magnetic  BL). 

Following  the  work  of  Godon  &  Sion  (2005),  we  use  the 
spin-down  time  rspin  =  116  s  and  the  rotation  (or  dynamical) 
time  rrot  =  Ircr/vKir)  =  25.3  s  to  derive  the  viscous  time 
in  the  BL,  r„  =  T2pin/Trot  =  532  s.  This  gives  a  BL  viscosity 
v  =  <5bl/t„  =  4.8  x  1013  cm2  s'1,  where  we  have  assumed 
Mw(j  =  1.29  Mq,  R,xl\  =  4000  Km,  and  a  BL  size  <5BL  given  by 
the  BL  radius  rB l  =  (1  +  <5BL)RW( i  =  1.4Rw(j  (Popham  1999). 
In  the  alpha  viscosity  prescription  v  =  acsH,  the  value  of  the 
alpha  parameter  is  then  simply  a  =  v/(csH )  ~  0.003,  where 
we  assumed  cs  &  108  cm  s'1  (for  T  ~  108  K)  and  the  vertical 
thickness  of  the  BL  H  ~  0.4/\\v(|  (Popham  1999).  This  value 
of  a  for  the  BL  of  RU  Peg  and  the  one  derived  for  the  BL  of 
VW  Hyi  (0.009;  Godon  &  Sion  2005)  indicate  that  the  alpha 


viscosity  parameter  in  the  BL  (q:Bl)  is  much  smaller  than  in  the 
disk  (adisk;  at  least  in  quiescent  DNe).  This  is  consistent  with  the 
analytical  estimates  of  Shakura  &  Sunyaev  (1988)  and  Godon 
(1995)  which  obtained  a^L  adisk-  Since  the  source  of  the 
BL  viscosity  is  unknown  (Inogamov  &  Sunyaev  1999;  Popham 
&  Sunyaev  2001),  this  result  is  important  for  future  theoretical 
work  investigating  the  source  of  the  viscosity  in  the  BL. 

P.G.  wishes  to  thank  Bill  Blair  for  his  kind  hospitality  at  the 
Johns  Hopkins  University  where  part  of  this  work  was  carried 
out.  Except  for  the  symbol  and  temperature  mark  of  RU  Peg, 
Figure  6  was  taken  from  Pandel  et  al.  (2005),  who  kindly  agreed 
that  we  reproduce  Figure  4  from  their  original  work.  This 
work  is  based  on  observations  obtained  with  XMM-Newton, 
an  ESA  science  mission  with  instruments  and  contributions 
directly  funded  by  ESA  member  states  and  by  NASA.  Support 
for  this  work  was  provided  by  NASA  through  grant  numbers 
NNX08AX43G  Q XMM-Newton  A07)  to  Villanova  University. 
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